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The effects of 5 mesons on the dynamical instabilities of cold and warm nuclear and stellar matter 
at subsaturation densities are studied in the framework of relativistic mean-field hadron models 
(NL3, NLp and NLp5) with the inclusion of the electromagnetic field. The distillation effect and the 
spinodals for all the models considered are discussed. The crust-core transition density and pressure 
are obtained as a function of temperature for /3-equilibrium matter with and without neutrino 
trapping. An estimation of the size of the clusters formed in the non- homogeneous phase is made. 
It is shown that cluster sizes increase with temperature. The effects of the (5-meson on the instability 
region are larger for low temperatures, very asymmetric matter and densities close to the spinodal 
surface. It increases the distillation effect above ~ 0.4po and has the opposite effect below that 
density. 
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INTRODUCTION 



Compact stars are to date believed to be made of inner 
layers enclosed in a crust and possibly, a shallow atmo- 
sphere Q . The information we can collect from the inte- 
rior of the star has necessarily crossed its crust. There- 
fore, a better understanding of these compact objects de- 
mands a precise description of their envelope. The inner 
crust of compact stars is described by a liquid-gas phase 
transition of asymmetric nuclear matter in the presence 
of electrons. In this context, the isospin content of the 
liquid and the gas phases can play an important role in 
transport phenomena and as consequence, account for 
the energy losses from neutrino emissivity [3, Q • 
Several theoretical and experimental efforts have been set 
in order to more precisely describe the physics involved 
in these phenomena, and there is currently considerable 
research and a number of experiments being done in the 
search for a better description of nuclear matter under 
exotic conditions (namely of density, temperature and 
asymmetry) . 

The authors of Ref. have stressed the importance 

of including the scalar isovector virtual (5(ao(980)) field in 
hadronic effective field theories when asymmetric nuclear 
matter is studied. Its presence introduces in the isovector 
channel the structure of relativistic interactions, where a 
balance between a scalar (attractive) and a vector (re- 
pulsive) potential exists. The (5, and p mesons give rise 
to the corresponding attractive and repulsive potentials 
in the isovector channel. The introduction of 5 mesons 
affects the behavior of the system at high densities, when 
its contribution is reduced leading to a harder equation of 
state (EOS) and at subsaturation density when its con- 
tribution is larger leading to a softer symmetry energy. 
The effects of the inclusion of 5 mesons on the properties 
of compact stars, such as mass, radius and strangeness 
content were discussed in d, 0] ■ In Q the 5-meson effect 
on the density dependence of the symmetry energy, on 



the nucleon transport effects and on resonance and par- 
ticle production around the threshold was studied. At 
subsaturation densities, the effect of the 5 meson was in- 
vestigated in the pasta formation at the inner crust of 
a compact star [91], the extension of the spinodal region 
[lo| or the spinodal region in the presence of very strong 
magnetic fields as the ones that might occur in magnetars 



The main goal of this work is to study the effect of 
5 mesons on the dynamical instabilities and phase tran- 
sitions in nuclear matter within the framework of rela- 
tivistic models investigating the effects of 5 mesons on 
matter at finite temperatures and under the conditions 
of isospin asymmetry and density expected in the inner 
crust. It is important to try these models at different 
conditions in order to better understand the properties 
of the neutron star crust. Some works have treated clus- 
terization of matter both in the non-relativistic [13, 
and relativistic [ij, [l^ contexts. The relativistic Vlasov 
equation formalism applied in the last two papers will 
be used in the present work. In [Tsj warm nuclear mat- 
ter matter was studied within NL3 [l^. We use three 
different models - NL3 NLp^ and NLp [J - to 

try to understand the role of 5 mesons in that context. 
We will consider both neutral neutron-proton-electron 
{npe~) neutrino free matter in /3-equilibrium at zero tem- 
perature and npe~ matter with trapped neutrinos in /3- 
equilibrium for a lepton fraction Yl^ — 0.4. We will also 



stablish a comparison with the results discussed in [17 1 
for non-relativistic Skyrme models and the NL3 relativis- 
itc model. 



In Sec. [ni we show the formalism we use and in Sec. 
Hm we present and discuss some of the results obtained 
and finally, in Sec. IIVI some conclusions are taken. 
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II. THE VLASOV EQUATION FORMALISM 

We use the relativistic non-linear Walecka model 
(NLWM) in the Mean-Field Approximation, within the 
Vlasov formalism to study nuclear collective modes of 
npe~ matter at finite temperature [isj . 

We consider a system of baryons, with mass M inter- 
acting with and through an isoscalar-scalar field (j) with 
mass nis, an isoscalar- vector field with mass rriy, an 
isovector-scalar field A with mass ms and an isovector- 
vector field with mass nip. We also include a system of 
electrons with mass mg. Protons and electrons interact 
through the electromagnetic field A^. The Lagrangian 
density reads: 

i—p,n 

where the nuclcon Lagrangian reads 

= [7^zZ?^ - M*] , 

with 

iD*" = id" - - • b^ - eA'^i-t^ , 
^ 2 2 



M* 



M - gs 



9 ST 



and the electron Lagrangian is given by 

The isoscalar part is associated with the scalar sigma 
(cr) field (j>, and the vector omega (uj) field Vp, whereas 
the isospin dependence comes from the isovector-scalar 
delta {S) field A*, and the isovector- vector rho (p) field 
6^ (where fj, stands for the four dimensional space-time 
indices and i the three-dimensional isospin direction in- 
dex). The associated Lagrangians are: 



= 

c — 

Ca = 
where f2,,, 



\ 



2 ±2 



12^'^ 



1 
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f 
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9uK 



(1) 



gp{hf, X b^) and ^ df,A^ 

The model comprises the following parameters: four 
coupling constants gs, gv, gs and gp of the mesons to the 
nucleons, the bare nucleon mass M, the electron mass 
rrie, the masses of t he meson s, the electromagnetic cou- 
pling constant e — ■y/47r/f37 and the self-interacting cou- 
pling constants k and A. In this Lagrangian density, t is 
the isospin operator. 



Model 


pa 


E/A 


K 


Esyui 


L 


Ksym 




NL3 


0.148 


-16.240 


269.936 


37.344 


118.320 


100.525 


-696.129 


NLp 


0.160 


-16.051 


239.884 


30.335 


84.510 


3.328 


-340.293 


NLp(5 


0.160 


-16.051 


239.884 


30.711 


102.664 


127.246 


-290.187 



TABLE I: Nuclear matter properties at saturation density, 
po- All quantities are in MeV, except for po, given in fm~'^. 



We have used the set of constants I and II, identified 
here as NLp and NLpi5, respectively, taken from [3| and 
also the NL3 parametrization [l^ . In the first two cases, 
the saturation density that we refer as po is 0.160 fm~^. 
For NL3, po = 0.148 fm-^. 

Table U shows some nuclear matter properties at the 
saturation density for these models: the binding energy 
per nucleon, the incompressibility coefficient, the sym- 
metry energy, the symmetry energy slope, the symmetry 



energy curvature and Kt — Ksym — 



Qo 
K 



L. 



All these properties, except Kgym and K^, present 
lower values for NLp and NLpS, as compared to NL3. 
The two parametrizations, NLp and NLp(5, have L values 
that are within the expected range estimated from re- 
cent experimental constraints: L = 88±25 MeV [l9l.[20l|. 
However, the values of Kt fall out of the Kr = — 550±100 
MeV mi, HI range. 

We denote by 

/(r,P,^)± = diag(/p±, /„±, /e±) 

the distribution functions of particles (+) at position r, 
instant t and momentum p and of antiparticles (-) at 
position r, instant t and momentum — p, and by 



h± = diag(/ip±, /i„±, he±) 
the corresponding one-body haniiltonian, where 



For protons and neutrons, i = p,n, we have 



(2) 



(3) 



M* = M - gsci) - gsnAs 

with Ti = I (protons) or -1 (neutrons). For electrons, 
I = e, we have 



Ve = -eA, 



m: 



The time evolution of the distribution functions is de- 
scribed by the Vlasov equation 



at 



Pi n, e. 



(4) 



where {, } denote the Poisson brackets. 
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From Hamilton's equations we derive the equations de- 
scribing the time evolution of the fields (j), V^, A^, the 
third component of the p-field, 63 = (60, b), and the third 
isospin component, A3, of the A field, which are given 
in Appendix [Al 

The state which minimizes the energy of asymmetric 
nuclear matter is characterized by the distribution func- 
tions 



1 + e(<^o»T'^.)/T' 



with 

and 

with 



M* 



-9vVo-Y^i"o-eAo — Y 



TOe± 



1 _(_ g(£OeTA'e)/T 



and by the constant mesonic fields, which obey the fol- 
lowing equations: ml(j)'^'' + f'/>o'^''^ + f0o°^^ = 9spf\ 

(0) _ „ _2l(0) _ ffp -(0) _2a(0) _ „_(o) 



r(") _ ui") _ ,1(0) 



,2 AO) _ gp AO) 

K"'=b^'^ = A^' =Af =0. 

Collective modes in the present approach correspond 
to small oscillations around the equilibrium state. These 
small deviations are described by the linearized equations 
of motion and, therefore, collective modes are given as 
solutions of those equations. To construct them, let us 
define: 



fi± 
4>0 
A3 

Vo 
bo 



_ f(0) 
— JOi± ' 

= 

= v^^ 

= SAo, 



S<P, 

-SA3, 

5bo, 
A, = 



As in [I^, [23, [2J], we express the fluctuations of the 
distribution functions in terms of the generating func- 
tions: 



S±ir,p,t) = diag(S'p±, S'„±, Se±) , 



such that 



5f± = {S±Jo.±}^{S±,p'}^^°^ 



In terms of the generating functions, the linearized 
Vlasov equations for (5/i±, 



d6fi± 
dt 



{Sfi±, ho,±} + {foi±,6hi±} = 



are equivalent to the following time-evolution equations: 



i± 



dt 



+ {Si±, hQi±} 



Shi± = T 



dS, 



e± 



dt 



{Se±,hoe±} — 5he± — 



MoT 



p-(5A 



eoe 



where 



g^5Vo + Sbo + e SAq, 

;(0) 



with hoi± = ±eoi + V^^' and /ioe± = ±eoe- 

Of particular interest on account of their physical rele- 
vance are the longitudinal modes, with momentum k and 
frequency uj, described by the ansatz 



/ 5,±(r,p,t) \ 



V 



<5A3 



/ Si±ip,cose)\ 



J 



V 



<^A3. 

5ii 



^i((^t-k-r) 



where j = p, n, e, ^ = V, b, A represents the vector- 
meson fields and 9 is the angle between p and k. For 
these modes, we get SV^ = SV^ = , db^ = Sb-^ = and 

6AZ = SAy^o. 

Calling ,5A3c. = SA^, 6V^ = SV^, Sb'^ = Sb^ and 
Mf, = SA^, we will have dV,^z = §y^^QK'^t-k■r) ^nd 
(5Voi = 6V^^e^^'^*~^'^\ In the Appendix we present the 
equations for the fields (eqs. (|Aip - (|A8p ). the equations 
of motion for the fluctuations (eqs. (|B4p - (jB10[) ). the so- 
lutions for the eigenmodes (jCip . as well as defining the 
coefficients aij. There we also show the amplitudes p^i 
fea. (IC4p ) and p^i (eq. (jC3l) ') obtained from the dispersion 
relation (eq. IC2p . 



III. 



NUMERICAL RESULTS AND 
DISCUSSIONS 



From the dispersion relation (see Appendix [Cl eq. (|C2p ) 
we have obtained the dynamical spinodal surfaces, char- 
acterized by a zero frequency, for different temperatures 
and momentum transfer, the ratios of the proton to neu- 
tron transitions densities (see eg. ljCSp ) which allow us to 
discuss the distillation effect. We also show the unstable 
modes with the largest energy modulus, which defines 
the mode that drives the system to a non-homogeneous 
phase and, therefore, gives an estimation of the size of 
the clusters formed in the phase transition. In the se- 
quel, we will discuss the effect of 6 mesons in these three 
physical quantities. 
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Figure [T] shows the ratio of the proton to neutron den- 
sity fluctuations as a function of the transferred momen- 
tum at T = 5 and 10 MeV. At T = 5 MeV and for 
the largest density considered, i.e., p — 0.5 poi the 5 me- 
son gives rise to a larger distillation effect. However, for 
p = 0.3po the effect of S is negligible and for p = 0.2po 
it reduces the distillation effect. This trend is not modi- 
fied with temperature, except that the instability for the 
larger densities may disappear. 

This effect is better seen from Fig. [2] where the ra- 
tio of the proton to neutron density fluctuations are also 
plotted as a function of the density for several values of 
the transferred momentum. As in Fig. [U we see that 
for a density above (below) p ~ 0.06 fm^'', the S meson 
increases (decreases) the distillation effect. This will af- 
fect the constitution of the crust: in the region closer to 
the inner edge the clusters are more proton rich, whereas 
the gas is more neutron-rich than the prevision without 
considering S mesons. For lower densities, the opposite 
happens. 

Figure [3] shows the subsaturation instability regions at 
T = 7 and 10 MeV for NLpS, NLp and NL3. It can be 
seen that these regions become smaller with the increase 
of temperature. The presence of the S meson also reduces 
slightly the spinodal region for very isospin asymmetric 
matter. For NLpS and NLp, the envelope of the spinodal 
region is characterized by fc ~ 100 MeV. On the other 
hand, NL3 has a smaller unstable region as compared 
with the other two parametrizations and the spinodal 
envelope is obtained for a smaller k, k ^ 75 MeV. The 
different k dependencies originate in the finite-range part 
of the nuclear force as discussed in . 

Figure [4] shows the approximate envelope of the spin- 
odal regions for different temperatures and the three 
models considered. The lines that cross the neutron- 
rich part of the spinodal are the EoS for /3-equilibrium 
matter both without neutrinos and considering neutrino- 
trapping with a lepton fraction = 0.4, as indicated in 
the figure. The non-homogeneous phase describing the 
crust corresponds to the EoS inside the spinodal. We 
notice that NLp and NLpi5 are almost coincident and 
only differ at large isospin asymmetries. We can see 
that the EoS for /9-equilibrium without neutrinos cross 
the spinodal region at T = MeV, and small tempera- 
tures. At T > 3 MeV, the non-homogeneous phase dis- 
appears in neutrino-free stellar matter. It means that 
after the neutrino outflow, a non homogenous phase ap- 
pears at the crust of the star when the star's surface has 
cooled down to a temperature below ~ 3 MeV. For mat- 
ter with trapped neutrinos, the EoS cross the spinodal 
for T < 11.6 MeV, considering the NL3 parametrization 
and T < 13.2 MeV, for the NLpS and NLp parametriza- 
tions. The non- homogeneous phase in the crust of the 
proto-neutron star will, therefore, exist if the tempera- 
ture of the crust is not larger than ^ 12 MeV being this 
critically model dependent. 

In Fig. Othe transition densities at the crust-core tran- 
sition for /3-equilibrium stellar matter with trapped neu- 



trinos (Yl = 0.4) are plotted as a function of temper- 
ature. Those values, as well as the corresponding pres- 
sures for different temperatures are also shown in Table 
im The crust-core transition densities decrease with tem- 
perature. Also, the pressure shows a similar behaviour, 
except for NL3 where it slightly increases. Pressures are 
slightly larger in the presence of 6 mesons. 

In T7| , the authors compare predictions from both 
nonrelativistic (Skyrme forces) and Relativistic Mean- 
Field (RMF) and density-dependent models (DDM). For 
Skyrme models, the crossing density for matter with 
trapped neutrinos at T=0 occurs 0.088 fm-3, except 
for two special cases, just slightly larger than most of 
the relativistic models studied and slightly below NLp 
and NLp(5. The critical temperature above which the in- 
stabilities disappear is ^ 12 MeV in agreement with the 
temperature obtained for NLp and NLpS. 



T[MeV] 


Pcros 3 [fm 


P[MeV.fm-^] 




NL3 


NLp5 


NLp 


NL3 


NLp<5 


NLp 





0.081 


0.091 


0.092 


0.995 


1.254 


1.242 


2 


0.079 


0.092 


0.092 


0.969 


1.256 


1.245 


5 


0.078 


0.088 


0.089 


1.087 


1.250 


1.242 


7 


0.073 


0.084 


0.085 


1.105 


1.249 


1.244 


10 


0.060 


0.075 


0.076 


0.971 


1.244 


1.243 


12 




0.064 


0.065 




1.122 


1.124 


13 




0.053 


0.055 




0.980 


0.982 



TABLE II: The crossing densities and pressures at the crust- 
core transition for /3-equiIibrium stellar matter with trapped 
neutrinos {Yl = 0.4) for several temperature and the three 
models: NL3, NLp5 and NLp. 



The unstable modes of the system are calculated by 
replacing the frequency a; by zF, where F defines the ex- 
ponential growth rate of the instabilities. We take the 
mode with the largest growth rate as the most unstable 
amongst all and, therefore, the one that drives the system 
to the formation of instabilities. Its half-wavelenght de- 
fines the most probable size of the clusters (liquid) formed 
in the mixed (liquid-gas) phase In Fig. [5] we plot 

the estimated size of the clusters in /3-equilibrium mat- 
ter without trapped neutrinos at T=0 MeV for the three 
models considered. It is shown that NL3 predicts larger 
clusters as compared with the other two models, and for 
a smaller density range. In Fig. [7] we also plot the es- 
timated size of the clusters, with trapped neutrinos and 
for different temperatures: T=5, 10 and 12 MeV. When 
neutrinos are included and the temperature is increased, 
there is a decrease in the density range for clusterization 
and an increase in the size of the clusters. It is interest- 
ing to notice the effect of temperature on the size of the 
clusters which does not agree with the conclusions of Ref. 
[2^ which predicts a small reduction of the cluster size 
with temperature (see Fig 7 of [1^). There, a simple ap- 
proach was used to determine the clusters: a zero surface 
thickness ansatz for the clusters was used with a non-self 



5 



T 




k (MeV) k (MeV) 

FIG. 1; (Color online) The ratio of the proton over the neutron density fluctuations plotted for yp — 0.2 and p — 0.2po (black), 
p — 0.3po (green) and p = 0.5po (blue) at T = 5 (a) and T — 10 MeV (b) as a function of the transferred momentum for NLp(5 
(dot-dashed) and NLp (continuous line). 




FIG. 2: (Color online) The ratio of the proton over the neutron density fluctuations plotted for i/p = 0.2 and fc = 10 MeV 
(black), fc = 50 MeV (green) and k = 100 MeV (blue) at T = 5 (a) and T = 10 MeV (b) as a function of the density for NLp<5 
(dot-dashed) and NLp (continuous line). 



consistent surface energy. It is, therefore, important to 
understand this point by describing the clusterized phase 
within a finite temperature self-consistent Thomas- Fermi 
calculation. 

At T=12 MeV the NL3 parameter set predicts no clus- 
ters. Furthermore, this model shows larger clusters and 
a smaller density range of the non-homogeneous phase 
as compared with the other two. At the transition den- 
sities, NLpS and NLp present clusters with the double 
of the size as compared with the intermediate densities. 
For T=10 MeV, Skyrme models predict a smaller range 
of the instability re gion and also smaller clusters at the 
transition densities as compared to NLp5 and NLp. 
Globally, the effect of the temperature is reducing the 
instability region and increasing the cluster size. 



IV. CONCLUSIONS 

In the present work we have studied the dynami- 
cal instabilities and phase transitions in nuclear mat- 
ter within the framework of the relativistic non-linear 



Walecka model (NLWM) in the Mean-Field Approxima- 
tion, using the Vlasov formalism. Moreover, we have 
used three different parametrizations, NL3 flG"!, NLp5 
and NLp 4] , in order to understand the role of S mesons 
in the crust-core phase transition of a cold and a warm 
neutron star. 

Summarizing, the S meson has a larger effect in 
neutron-rich matter, at larger densities and lower tem- 
peratures. For densities p > O.OGfm""^, the S meson in- 
creases the distillation effect and the contrary occurs be- 
low that density. In addition to it, the distillation effect 
decreases with temperature. 

Besides the distillation effect we have also tested the 
behaviour of the spinodal regions in the presence of S 
mesons and with the increase of temperature. It was 
shown that the spinodal regions become smaller at tem- 
perature increase, and slightly smaller with the inclusion 
of the S meson. The model NL3 has a smaller envelope 
of the instabilities. A larger spinodal for NLpS and NLp 
is related to the finite range of the nuclear force in this 
model. 

In order to study the presence of a non-homogeneous 
phase in /3-equilibrium stellar matter we have determined 
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FIG. 3: (Color online) Spinodal curves for fc = 30 MeV (black), fc = 75 and 100 MeV (blue), k = 150 and 200 MeV (green) at 
T = 7 ((a) and (b)) and T = 10 MeV ((c) and (d)) for NLpJ (dot-dashed), NLp (continuous) and NL3 (dashed line). 
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FIG. 4: (Color online) Spinodals for fc = 75 (black) and fc = 100 MeV (green) for NL3 (dashed), NLp5 (dot-dashed) and NLp 
(continuous line) at T = (a), T = 10 (b) and T = 12 MeV (c). The lines that cross the spinodals are the EoS in /3-equilibrium 
for npe~ (neutrino-free) and with trapped neutrinos for a lepton fraction Yl = 0.4 (red and blue, respectively), as indicated. 



the crossing density of the spinodals with the EoS for 
two cases: with neutrinos, where the constant leptonic 
fraction was Yl = 0.4 and in npe~ matter. For the first 
case, the EoS crosses the spinodal for T < 11.6 MeV, 
considering the NL3 parametrization and T < 13.2 MeV, 
for the NLp6 and NLp parametrizations. For the second 
case, the EoS crosses the spinodal region only for T < 3 



MeV. Around T > 3 MeV, the non homogeneous phase 
disappears. Also, the unstable regions are larger for NLpJ 
and NLp, compared with NL3 and it was shown that the 
6 meson slightly decreases the instabilities. 

A quantity of concern that has recently been given 
more attention is the core-crust transition density. In 
[2^ the authors relate this density with the fraction of 



7 
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FIG. 5: (Color online) The crossing densities at the crust- 
core transition for /3-equilibrium stellar matter with trapped 
neutrinos {Yl = 0.4) for NL3 (blue, dashed), NLpS (green, 
dot-dashed) and NLp (black, continuous line) as a function of 
temperature. 
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FIG. 6: (Color online) Size of the clusters in /^-equilibrium in 
neutrino-free matter at T = MeV for NL3 (blue, dashed), 
NLp(5 (green, dot-dashed) and NLp (black, continuous line). 



moment of inertia of compact stars. In [21] a compila- 
tion of neutron skin thicknesses of several nuclei brings a 
relation between the symmetry energy coefficients of the 
EoS of nuclear matter and of nuclei described by several 
different models, what allows to estimate the core-crust 
transition density, which is reported to be ~ 0.095 fm~'^. 
The values obtained in the present work at zero tempera- 
ture are slightly lower than this value and extremely close 
to the findings in [fl] where a Thomas Fermi calculation 
for the pasta phase was carried out. 

We have determined the crust-core phase transition 
densities and the corresponding pressures as a function 
of temperature. These two quantities decrease with the 
temperature for NLp and NLp(5 above T=2 MeV. For 
NL3, the pressure slightly increases for T< 7 MeV and 
decreases above that temperature. The pressures are 
slightly larger in the presence of S mesons. The inclu- 
sion of i5 only slightly reduces the crossing densities and 
increases the pressure at transition which is the main 
quantity defining the fraction of moment of inertia cor- 
responding to the crust [26j . 



The estimated size of the clusters was determined as a 
function of temperature. NL3 shows larger clusters and 
a smaller range of the non-homogeneous phase as com- 
pared with the other two. At the transition densities, 
NLp and NLpS predict clusters with double size as com- 
pared with an intermediate density. It was also shown 
that the size of the clusters increases with temperature. 
The main effect of S mesons on the size of the clusters oc- 
curs precisely at the crust-core transition: it reduces the 
transition density and the steep increase into larger clus- 
ters occurs at lower densities. A careful determination of 
the size of the clusters for the temperature at which neu- 
trino trapping occurs is important to understand whether 
this may affect the interaction of the neutrinos with nu- 
clear matter and favour the reactivation of the explosion 
in the dynamics of a supernova. 

It has been shown that at finite temperature the be- 
haviour of very asymmetric warm nuclear matter depends 
on the model, and that 5 mesons have a stronger effect 
at the crust-core transition density. This suggests that 
better constraints at finite temperature on the EoS are 
required. 



APPENDIX A: EQUATIONS FOR THE FIELDS 



- V^Ag + m^Ag =55P3.(r,t) , (A2) 

-V^Vo+ mlVo = gwo(r, t) , (A3) 
-V^V^ + mlV^ ^ g.Mv.t) , (A4) 



V%o + mlbo ^ ^j3,o{r,t) , (A5) 
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V^fe^-^m^fe^ = ^J3,.(r,t), (A6) 
V^Ao = e[jop(r,0 -ioe(r,0] , 



V^A, = efc(r,i) - j,e(r,t)] , 



(A7) 
(A8) 



where the scalar density is 



i—p^n 
= Psp + Ps 



(27r)3 ^-^"^ ^' ^ 
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(C) T=12 MeV 
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FIG. 7: (Color online) Size of the clusters in /3-equilibrium with neutrino-trapping (Yl — 0.4) at T = 5 (a), T = 10 (b) and 
T = 12 MeV (c) for NL3 (blue, dashed), 'NLpS (green, dot-dashed) and NLp (black, continuous line). 



and the isovector density is 



APPENDIX B: EQUATIONS OF MOTION 



P3.(r,i) = 2 ^ / T^r,(/,+ (r,p,i) + /,_(r,p,i))^ 

i—p,n 
Psp Psn- 

The components of the baryonic four-current density are 
jo(r,t) = J (|;^(/'+(r,P,t)-/.-(r,P,t)) 



Pp+ Pn , 



j(r,t) = 2 ^ I ^(/,(r,p,t) + /,_(r,p,i)) 



joe(r,i) = 2 J ^^-p- (/e+(r,p,t) - /e_(r,p,t)) , 



Ur,t)=2 J ^0j(/,+ (r,p,t) + /e_(r,p,t) 



p + eA 



and the components of the isovector four-current density 
are 

i—p^n 
~ Pp Pn : 



_ (27r)3 e, 

i—p,n 

X T, (/,+ (r,p,i) + /,_(r,p,i)) , 



with = -y/(p - Vif + Mf ,i^p,n and 
ie = v/(P + eA)2 +m2 . 



From the continuity equation for the density currents, 
we get for the components of the vector fields 

LuSbl = kSb^, 
uj6A" = kSA^. 



(Bl) 
(B2) 
(B3) 



Defining 

Spsi = 



5pi = 



with 



cfip pcos0 
(27r)3 eo, 

dfai± 1 



SI 



J, 



dfoi+ 


+ si_ 


dfoi- 


dp"^ 


dp"^ 


dfoi+ 




dfoi- 


dp'^ 




dp"^ 



dp^ 2T€a 
the equations of motion read 
kp cos 



/0'i±(/0i± - 1) 



kp COS t 



gs Ti M, 



5! 



It 



uj p COS 9 
k eoe 



5Al, (B4) 



9sM* 



- fc 



It 



w p COS 6* 



5Vl\ (B5) 



^ ml „ff) 54>uj = 4:igsk ^ M*6ps 

i—p,n 



gs dpl^ 6 A 



2 



;2-fc2 



i|gyy)(5A^ = Aigsk ^ TiM*Sps 

i—p,n 

+ gsgsdpl^Scj)^, 



t—Pj7l 



;2 - fc2 _ ™2 



) t^^" = 2igpk ^ r,5p,, 
_ k^) SAl = Aiek J2 



(B6) 

(B7) 
(B8) 

(B9) 

(BIO) 



2— p,e 
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where 



and 



A 



ml + K(j)o + '-4>f^ + gldp° 



(27r)3 ^ 



d?pn {foi+ + foi-) 



1 Af*2 



(27r)2rw^ 



(27r) 



n^T \kJ ' 



p/e 



dx- 



p/t 



= ± 



2 — V ijj In 



uj-p/e 



e'%^(e)/o^T(/o.T-l)*. 



A/* 



APPENDIX C: SOLUTIONS FOR THE 
EIGENMODES AND THE DISPERSION 
RELATION 



The sohitions of Eqs. (|Bl[) - (|B10p form a complete set 
of eigenmodes that may be used to construct a general 
solution for an arbitrary longitudinal pertubation. Sub- 
stituting the set of equations (|B6P - (|B10p into (|B4p and 
()B5[) we get a set of equations for the unknowns 5^^, 
which lead to the following matrix equation: 



(CI) 



The amplitudes pf^ (eq. (|C4l) ) and p^^ (eq.jGSj) are 
functions of the quantities Sl^j^. The dispersion relation 
is obtained from the determinant of the matrix of the 
coefficients, namely, 



^ ail 
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Puje 
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Det{a,j) = 0. 
We define the following quantities 



(C2) 



pcos9 



Grf,, — 



9sM* 



— 



[(02-C^2)(^2_^2) 



9s95dp°3s 
fc2 



fc2 



27r2T 

+ Gs, {{u,^ ~ cj^,)Gs^ + D,G^^)] 



(-j {1-u; ), = ^{k 



2 I 2\ 



rOO rp/e 

^Zi,± = / <^"d,e / dxxSl±{x,p) fo,^± (/oi.± - 1) ; 

JM* J-ph 



Pu>i 

- 4" 



I = p,n,e 



I = p,n . 



(C3) 
(C4) 



The coefhcients defined as: 



On 
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flu 
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jOp _ Mp 



oi2 = cr( 

-(c„ - Cp) [iX + I'J-) , 015 = + I'J_ 

c7 {C+ " C-) , 022 = 1 + cT (/*; - C-) , 

-(C„+Cp)(/1';+/I'i), 025=0, 

+C7 + I'J-) , 032 = + I'J-) , 

l-(c. + Cp + c,)(e+-/^^„), 

-(C« - Cp) (/^^ - ) , 035 = Ce 

np / rln I Tin ^ „ _ / rln , T^ln \ 



2p _ x2p 



2n 7-2n ' 



1 - (C. + Cp) - ) , 

052 = 054 = 0, 

Ce (e; - C-) , 055 = 1 - Ce (e+ - 

The ratios of the amplitudes are given by: 
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Pun 

with 

Opp 



; 1 1 1 1 1 



OllOp 



ai2a„ 



Oi4 



oiiflpp + ai2a„p + ai3 — 015053/055 



(C5) 



022043 — 023042 
021O42 — 022041 

044021 — 041024 
041022 — 042021 ' 

Pue 
Pup 



O22O44 — O24O42 
O21O42 — O22O41 

O43O21 — O41O23 
O41O22 — O42O21 



O53 
O55 



(C6) 



10 



ACKNOWLEDGMENTS SFRH/BPD/29057/2006 and CERN/FP/83505/2008. 

This work was partially supported by FCT 
(Portugal) under grants PTDC/FIS/64707/2006, 



[1] M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis, J. M. 

Lattimer, and R. Knorrcri, Phys. Rep. 280, 1 (1997). 
[2] R. F. Sawyer, Phys. Rev. D 11, 2740 (1975). 
[3] N. Iwamoto and C. J. Pethick, Phys. Rev. D 25, 313 

(1982). 

[4] B. Liu, V. Greco, V. Baran, M. Colonna, and M. Di Toro, 

Phys. Rev. C 65, 045201 (2002). 
[5] V. Greco, M. Colonna, M. Di Toro, and F. Matera, Phys. 

Rev. C 67, 015203 (2003). 
[6] D. P. Menezes and G. Providencia, Phys. Rev. G 70, 

058801 (2004). 

[7] B. Liu, H. Guo, M. Di Toro, and V. Greco, Eur. Phys. J. 

25, 293 (2005). 
[8] T. Gaitanos, M. Di Toro, S. Typel, V. Baran, C. Fuchs, 
V. Greco, and H. H. Wolter, Nucl. Phys. A 732, 24 
(2004). 

[9] S. S. Avancini, L. Brito, J. R. MarineUi, D. P. Menezes, 

M. M. W. de Moraes, G. Providencia, and A. M. Santos, 

Phys. Rev. G 79, 035804 (2009). 
[10] S. S. Avancini, L. Brito, P. Chornaz, D. P. Menezes, and 

C. Providencia, Phys. Rev. C 74, 024317 (2006). 
[11] A. Rabhi, C. Providencia, and J. D. Providencia, 

arXiv:hep-th/0904.3261vl (2009). 
[12] P. Ghomaz, M. Golonna, and J. Randrup, Phys. Rep. 

389, 263 (2004). 
[13] V. M. Kolornietz and S. Shlomo, Phys. Rep. 390, 133 

(2004). 

[14] C. Providencia, L. Brito, S. S. Avancini, D. P. Menezes, 

and P. Ghomaz, Phys. Rev. G 73, 025805 (2006). 
[15] L. Brito, G. Providencia, A. M. Santos, S. S. Avancini, 



D. P. Menezes, and P. Ghomaz, Phys. Rev. G 74, 045801 

(2006). 

[16] G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. G 55, 
540 (1997). 

[17] G. Ducoin, G. Providencia, A. M. Santos, L. Brito, and 

P. Ghomaz, Phys. Rev. G 78, 055801 (2008). 
[18] M. Nielsen, G. Providencia, and J. da Providencia, Phys. 

Rev. G 44, 209 (1991). 
[19] J. Xu, L. W. Ghen, B. A. Li, and H. R. Ma, Phys. Rev. 

G 77, 014302 (2008). 
[20] J. Xu, L. W. Ghen, B. A. Li, and H. R. Ma, Phys. Rev. 

G 79, 035802 (2009). 
[21] T. Li, U. Garg, Y. Liu, R. Marks, B. K. Nayak, P. V. M. 

Rao, M. Fujiwara, H. Hashimoto, K. Kawase, K. Nakan- 

ishi, et al., Phys. Rev. Lett. 99, 162503 (2007). 
[22] U. Garg, T. Li, S. Okumura, H. Akimune, M. Fujiwara, 

M. Harakeh, H. Hashimoto, M. Itoh, Y. Iwao, T. Kawa- 

bata, et al., Nucl. Phys. A 788, 36 (2007). 
[23] S. Avancini, L. Brito, D. P. Menezes, and G. Providencia, 

Phys. Rev. G 71, 044323 (2005). 
[24] S. S. Avancini, L. Brito, D. P. Menezes, and 

G. Providencia, Phys. Rev. G 70, 015203 (2004). 
[25] S. S. Avancini, D. P. Menezes, M. D. Alloy, J. R. 

MarineUi, M. M. W. Moraes, and G. Providencia, Phys. 

Rev. G 78, 015802 (2008). 
[26] B. Link, R. I. Epstein, and J. M. Lattimer, Phys. Rev. 

Lett. 83, 3362 (1999). 
[27] M. Gentelles, X. Roca-Maza, X. Vifias, and M. Warda, 

Phys. Rev. Lett. 102, 122502 (2009). 



